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We use the newly released 182 Type la supernova data combined with the third-year Wilkinson 
Microwave Anisotropic Probe data (WMAP3) and large scale structure (LSS) information including 
SDSS and 2dFGRS to constrain the dark energy equation of state (EoS) as well as the curvature 
of universe Qk- Using the full dataset of Cosmic Microwave Background (CMB) and LSS rather 
than the shift parameter and linear growth factor, we make a Markov Chain Monte Carlo (MCMC) 
global fit, while paying particular attention to the dark energy perturbation. Parameterizing the 
O , EoS as WDE^a) = wo+w>i(l — a), we find the best fit of (wo, wi) is (—1.053, 0.944) and for wde(ci) = 

wo + mi sin(|7r ln(a)), the best fit for (wo,wi) is ( — 1.614, —1.046). We find that a fiat universe is 
a good approximation, namely, |fijc| > 0.06 has been excluded by 2a yet the inclusion of Qk can 
affect the measurement of DE parameters owing to their correlation and the present systematic 
effects of SNIa measurements. 
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I. INTRODUCTION 
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Dark energy (DE), the very power to drive universe's acceleration, is one of the most important issues in modern 
cosmology. Its existence was firstly revealed by the measurement of the relationship between redshift z and luminosity 
distance dj, of Type la supernova (SN Ia)[lJ. Dark Energy encodes its mystery in its equation of state (EoS) defined 
as the ratio of pressure over energy density thus DE models can be classified in terms of EoS@. 

The simplest candidate of dark energy is the cosmological constant (CC) whose EoS remains —1. Favored by 
current astronomical observations as it is, CC suffers from severe theoretical drawbacks such as the fine-tunning 
and coincidence problemQ. Alternative DE models with rolling scalar field, such as quintessence phantomQ, 
k-essence@, etc have been studied. The EoS of these models varies with cosmic time either above —1 or below — 1 
during evolution but the statement of "No-Go" theorem forbids it to cross the —1 boundary Q. Models where gravity 
is modified can also give these observed effects. 

Given our ignorance of the nature of dark energy, constraining the evolution of DE the EoS by cosmological 
observations is of great significance. Various methods have been used to constrain DE including parametric fitting^, 
H^UIli non-parametric reconstruction [ll|, etc. Interes ting ly, there exists some hint that the EoS of DE has crossed over 
k> \ —1 at least once from current observations [9l [Til. [l2l fl3f. which greatly challenges the above mentioned dark en ergy 
models, albeit the evidence is still marred by systematic effects. Quintom, whose EoS can smoothly cross — 
has attracted a lot of attention in the literature since its invention [l^. There have been many efforts in quintom 
model building, for example, double-scalar-field realization [lil HH] . a single scalar field with high derivative 17|, vector 
fields and so forth. 

A special and interesting example of quintom is Oscillating Quintom, whose EoS oscillates with time and crosses 
— 1 many times. The oscillating behavior in the EoS leads to oscillations in the Hubble constant and a recurrent 
universe. Oscillating Quintom is physically well motivated, since this scenario, to some extent, unifies early inflation 
and the current acceleration of the universe [l9j. In Ref . (2(j| . we have presented some preliminary studies on oscillating 
quintom. 

The nature of dark energy is a dominant factor of the fate of our universe. Another critical point is the curvature, 
fix- neutrino masses, which probably exists, albeit small, and the curvature, SIk- We concentrate, in this paper, on 
the correlation between VLk and the dark energy parameters. 1 

We use the newly released 182 SN la "Gold sample" (SN182)[H combined with WMAP3[H 2 and LSS information 
to constrain the evolution of the DE EoS and the curvature of universe. For DE EoS, we choose two parameterizations 



1 Other cosmological parameters also affect the probing of dark energy, such as neutrino mass [Sill and inflationary parameters |22 

2 Available at http://lambda.gsfc.nasa.gov/product/map/current/ 
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as in (JTJ) and ([2]) and we will address our motivation for such a choice in the next section. In our study, we treat the 
curvature of the universe £Ik as a free parameter rather than simply assuming a flat universe and make a Markov 
Chain Monte Carlo (MCMC) global fit based on Bayesian statistics. Paying particular attention to the dark energy 
perturbation especially when EoS crosses — lQ, we find the latest observations mildly favor quintom model however 
ACDM remains a good fit. We have also found that the inclusion of VIk can affect the determination of DE parameters 
significantly due to their correlation. 

We structure this paper as follows: after this introductory part, we propose our method and define the data set 
used in section II. The in section III we present our results and end up with discussion and comments. 



II. METHOD AND DATA 



To study the dynamical behavior of dark energy, we choose two kinds of parametrization of dark energy equation 
of state: 
I) 

wde{o) — w + wi(l - a) (1) 

II) 

WDE(a) — w + wi sin(w 2 hi(a)) (2) 

where a is the scale factor, wo denotes the EoS at present epoch and w± and W2 characterize the time evolution of 
DE. Parametrization I) is the most popular in literature since W\ simply equals to —dwoE{o) / da, which is the time 
derivative of «;nw(a)[2q]. Thus it is straightforward to study the dynamical behavior of DE. The physical motivation 
of parametrization II) is oscillating quintom. From @ we can see at low redshift, II) takes a form similar to I). At 
medium and high redshift, the EoS keeps oscillating. 



From the latest SN la paper[24(, one can find some hint of oscillating behavior of the EoS in their FIG. 10 where 
they use a quartic polynomial fit. Our sine function has the advantage of preserving the oscillating feature of the EoS 
at high redshift measured by the CMB data. For simplicity and focus on the study at lower redshift, we set W2 to be 
§7r in order to allow the EoS to evolve more than one period within the redshift range of to 2 where SN data are 
most robust. 

When using the MCMC global fitting strategy to constrain cosmological parameters, it is crucial to include dark 
energy perturbation. This issue has been realized by many researchers including the WMAP group [3, H, [2f| |27|. How- 
ever one cannot handle the dark energy perturbation when the parameterized EoS crosses —1 based on quintessence, 
phantom, k-essence and other non-crossing models. By virtue of quintom, the perturbation at the crossing points is 
continuous, thus we have proposed a sophisticated technique to treat dark energy perturbation in the whole parameter 
space, say, EoS> —1, < —1 and at the crossing pivots. For details of this method, we refer the readers to our previous 
companion paper @, [|[ . 

In this study, we have modified the publicly available Markov Chain Monte Carlo package C AMB / CosmoMC [28j 
to include the dark energy perturbation when the equation of state crosses — 1. 

The dark energy EoS and curvature of the universe fix can affect the determination of the geometry of our universe 
thus DE parameters are correlated with fix- Therefore, in our study, we relax the curvature of universe VLk as a free 
parameter rather than simply assuming a flat universe. We assume purely adiabatic initial conditions and set our 
most general parameter space as: 

P = cj c , 6 S , r, n K , w , w u n„ ln(10 10 A s )) (3) 

where lo^ = Vt^h 2 and lo c = VL c h 2 are the physical baryon and Cold Dark Matter densities relative to the critical 
density, Q s is the ratio (multiplied by 100) of the sound horizon to the angular diameter distance at decoupling, r 
is the optical depth to re-ionization, Qk = 1 — $l m — ^de is the spatial curvature, Wo,Wi portray the dynamical 
feature of dark energy as illustrated in ([1]) and ^j. A s and n s characterize the primordial scalar power spectrum. 
For the pivot of the primordial spectrum we set fc s o = 0.05Mpc _1 . Furthermore, we make use of the Hubble Space 
Telescope (HST) measurement of the Hubble parameter Hq = lOOh km s _1 Mpc _1 [2!| by multiplying the likelihood 
by a Gaussian likelihood function centered around h = 0.72 and with a standard deviation a = 0.08. We also impose 
a weak Gaussian prior on the baryon density f2f,/i 2 = 0.022 ± 0.002 (1 a) from Big Bang Nucleosynthesis [30]. The 
weak priors we take are as follows: r < 0.8, 0.5 < n s < 1.5, —0.3 < £Ik < 0.3, — 3 < wq < 3, — 5 < w\ < 5 3 and a 



We set the prior of uiq and wi broad enough to ensure the EoS can evolve in the whole parameter space. 
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TABLE 1. Constraints of dark energy equation of state and some background parameters when relaxing Qk as a free 
parameter (Left panel) and assuming a fiat universe (Right). For each case, we consider two forms of parametrization of dark 
energy EoS: Linear (w(a) = wo + Wi(l — a)) and Oscillating {w(a) — wo + »i sin(iU2 ln(o)), set W2 = §7r, see text for 
explanation). Best fit models, which give the minimum x 2 , an d the marginalized 2a errors are shown. All these constraints 

are from data combination of WMAP3 + SN182 + SDSS + 2dFGRS. 
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cosmic age tophat prior 10 Gyr < to < 20 Gyr. 

In our calculations, we have taken the total likelihood to be the products of the separate likelihoods £ of CMB, 
LSS and SNIa. In other words, defining x 2 = ^21og£, we get 

2 _ 2 ,2,2 { A\ 

Xtotal — XCMB T XhSS "+" XsNIa ■ V 4 J 

If the likelihood function is exactly Gaussian, x 2 coincides with the usual definition of % 2 up to an additive constant 
corresponding to the logarithm of the normalization factor of C. In the calculation of the likelihood from SNIa we have 
marginalized over the nuisance parameter [3l| . The supernova data we use are the "gold" set of 182 SNIa recently 
published by Riess et al in Ref.[24]. In the computation of CMB we have used the full dataset of the WMAP3 data 
with the routine for computing the likelihood supplied by the WMAP team [25j . For LSS information, we have used 
the 3D power spectrum of galaxies from the SDSS [IH and 2dFGRS[33f. To be conservative but more robust, in the 
fittings to the 3D power spectrum of galaxies from the SDSS, we have used the first 14 bins only, which are supposed 
to be well within the linear regime [34j . 

For each regular calculation, we run six independent chains comprising 150, 000 — 300, 000 chain elements and spend 
thousands of CPU hours to calculate on a cluster. The average acceptance rate is about 40%. We discard the first 
30% chain elements to be the "burn- in" process, test the convergence of the chains by Gelman and Rubin criteria [35j 
and find R — 1 of order 0.01, which is more conservative than the recommended value R — 1 < 0.1. 



III. RESULTS 



We summarize our main results in Table I. For all the combined data (SN182+WMAP3+SDSS+2dFGRS), we find 
that the flat universe is a good fit since the best fit value of VIk is —0.015 and —0.012 for DE parametrization I) and 
II) respectively and \£Ik\ > 0.06 has been excluded for more than 2a for both DE parameterizations. This can be seen 
graphically in FIG.l. From these 2-D contour plots of energy density of dark energy and matter from different data 
combination and different DE parameterizations, we find that the data of supernova-only favor a non-flat universe, 
however when CMB and LSS data are combined, a flat universe is preferred. 

For parametrization I), we find the best fit value of (wq, wi) to be (—1.149, 1.107) for a flat universe. When Q,k 
is freely relaxed, the best fit value of (wq,wi) is changed to (—1.053,0.944) and the error bars of nearly all the 
cosmological parameters have been enlarged. We find dark energy models whose EoS can cross —1 are mildly favored. 
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The 1-D posterior distribution of wo, Wi, £Ik and their 2-D correlation are shown in FIG. 2. In the wq — Q.k and 
W\ — fix panel, we find interesting correlation among curvature and DE parameters. This is expected since £Ik can 
contribute to luminosity distance via: 



d L (*) - 



1 + z 



=smn 



dz' 
E(z') 



(5) 



E(z) = ^ = [O ro (l + zf + n DE exp U f 1 + W[z ' i ) dz ; ) + n K (l + z) 2 ] 1/2 . (6) 
Uo \ Jo 1 + z J 

where sinn(y^|fc|ir)/-v/|fc[ = sin(a;), x, sinh(a;) if k = 1, 0, — 1. Furthermore, £Ik can modify the angular diameter 
distance to last scattering surface and the transfer function, which leaves imprints on the CMB and matter power 
spectrum. 

In the wq — wi panel of FIG. 2, the parameter space has been divided into four parts. The upper right and lower 
left parts denote for w > — 1 and w < — 1, the regions for quintessence and phantom models respectively. The other 
two parts dubbed "Quintom A(B)" represent models whose EoS can cross -1 during evolution. Quintom A crosses —1 
from upside down while Quintom B transits —1 from the opposite direction. We also plot the results when assuming 
a flat universe for comparison. We see the best fit model is in the region of Quintom A and the ACDM(the intersect 
of two dot dashed lines) is still a good fit. Moreover, relaxing £Ik enlarges the wo — w\ contour as expected. 

For parametrization II), again we find a small absolute value of Qk- The best fit models can cross -1 twice in the 
evolution. In FIG. 3, we show the correlation between and dark energy parameters, and we find the Quintom 
model whose EoS crosses -1 during evolution is preferred. 

We can see the dynamics of dark energy more clearly from FIG. 4. We show the best fit model and 2a errors of 
w(z) for the case of flat universe and relaxing Qk as a free parameter for the two DE parameterizations. The best fit 
model of each case is quintom-likc. 



IV. SUMMARY AND DISCUSSIONS 



In this paper we investigate the dynamics of dark energy and curvature of universe from the data of newly released 
182 supernova data combined with CMB and LSS information. Rather than assuming a flat universe, we relax Qk 
as a free parameter and make a MCMC global fit to measure the dark energy parameters as well as the curvature 
of universe. We find the model whose EoS can cross over —1 is favored for two dark energy parameterizations we 
considered in this work, albeit the ACDM model remains a good fit. A flat universe is preferred, namely, the best fit 
value of \Qk\ is smaller than 0.015 for the two DE parameterizations and \SIk\ > 0.06 has been excluded by more than 
2a. However, the correlation among dark energy parameters and Qk might not be neglected. Freeing fix enlarges 
the contours and even modifies the best fit value of dark energy parameters. For example, for parametrization II), 
relaxing tlx has changed the best fit value of (wo, Wi) from (—1.149, —0.525) to (—1.614, —1.046). This is because 
adding CIr can reduce the total \ 2 by 1.2 thus the global minimum moves to a deeper point. 

For CMB information, we use the full dataset of WMAP3, rather than the shift parameter. The shift parameter has 
the advantage of easy implcmentaton and much shorter calculation time, but it is not really "model independent" . It 
is derived from a fiducial ACDM model, thus it might lead to biased results, if one fits to models departing significantly 
from the fiducial model. Another drawback is that the shift parameter can merely offer part of the CMB information 
related to background parameters, thus it cannot constrain the pcrturbated dark energy. 

We will have a deeper and deeper understanding of dark energy with the accumulation of high quality cosmological 
observation especially for supernova data, such as future SNAP, ESSENCE, etc. To be bias-free as much as possible, 
it is better to add Hk into the global analysis, use full datasets of CMB and carefully treat dark energy perturbation. 
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FIG. 1: 2-D contour plots of energy density of matter and dark energy using latest astronomical observations. Grey: SN182; 
Blue: WMAP3; Red: WMAP3 + SN182 + SDSS + 2dFGRS. The dark and light shaded regions stand for 68% and 95% 
C.L. respectively. Different dynamical behaviors of dark energy have been considered. Left: vacuum energy (w — —1 forever); 
Middle: linearly growing (w(a) = Wq + Wi(l — a)); Right: Oscillating {w(a) — wq + W\ sin(u>2 ln(a))), with W2 = |tt. 
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FIG. 2: Constraints on the curvature Qk and dark energy parameters wo and wi when parameterizing DE EoS as w(a) = 
wo + wi(l — a). 1-D plots show the posterior distribution of £Ik,wo and i«i while the 2-D contour plots illustrate their correlation. 
The dark and light shaded area stand for 68% and 95% C.L. respectively. Blue dashed curves denote the case of flat universe 
for comparison. Different dark energy models can be distinguished from the wo — wi panel (See text). All these constraints are 
from data of WMAP3 + SN182 + SDSS + 2dFGRS. 



8 




FIG. 3: The same graphic convention as in FIG. 2 except for the parametrization of the EoS of DE: w(a) = Wo + wi sin(iU2 ln(a)). 
For simplicity, we have chosen W2 to be |7r. Note the variations of the error on for negative £1%, as a function of Wi, which 
reflect the oscillating feature of this parametrization of the EoS. 
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FIG. 4: Constraints of the time evolving of equation of state of DE using WMAP3 + SN182 + SDSS + 2dFGRS. Upper panel: 
w(a) — wo + wi(l — a); Lower panel: w(a) = U)q + wi sin(iU2 ln(a)). The cases of flat universe and treating Qk as a free 
parameter are both considered as illustrated in the plot. The solid blue lines denote the best fit models while the shaded areas 
illustrate the 2a errors. Red dashed lines show the cosmological constant boundary. 



